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We present in this article a novel Lagrangian measurement technique: an instrumented particle 
which continuously transmits the force/acceleration acting on it as it is advected in a flow. We 
develop signal processing methods to extract information on the flow from the acceleration signal 
transmitted by the particle. Notably, we are able to characterize the force acting on the particle 
and to identify the presence of a permanent large-scale vortex structure. Our technique provides a 
fast, robust and efficient tool to characterize flows, and it is particularly suited to obtain Lagrangian 
statistics along long trajectories or in cases where optical measurement techniques are not or hardly 
applicable. 

Turbulence is omnipresent in nature and in industry, and has received much attention for years. In 
the specific field of experimental fluid dynamics research, very significant progress has been achieved 
during the last decade with the advent of space and time resolved optical techniques based on high speed 
imaging [ ]. However, a direct resolution of the Eulerian flow pattern is still not always possible nor 
simple to carry out. In this context, Lagrangian techniques, in which the fields are monitored along the 
trajectories of particles, provide an interesting alternative [2, 3] with information about the small scales of 
turbulence (especially isotropy) and a huge focus on the particle's Lagrangian acceleration that directly 
reflects the turbulent forces exerted on the particles [4-9]. 

From an experimental point of view several inconveniences arise. In the Lagrangian frame one would 
like to collect long trajectories. However, even in confined flows it is difficult to track even just a few 
particles over a long time using the existing methods. For instance to use optical methods, the flow 
must be entirely observed and continuously recorded, something which is not yet possible. Apart from its 
implication for computing converged statistical quantities, several theories such as the fluctuation theorem 
necessitate long trajectories instead of many short ones. Another issue is the possible rotation of large 
particles in a flow, and the influence of this possible rotation on the dynamics of the particle. An optical 
technique following simultaneously particle position and absolute orientation in time has been recently 
developed [ ]. It shows in particular that for increasing turbulence, solid particles experience stronger 
rotation [11, 12]. The technique used in those experiments is not straightforward and needs careful 
calibration and synchronization, as well as an expensive set-up (high-speed cameras, strong illumination, 
etc) as well as time-consuming post-processing. Other common Lagrangian techniques, e.g. particle 
tracking velocimetry (PTV), generally do not allow a direct measurement of the possible rotation of the 
particle simultaneously with its translation. 

The experimental technique presented here was designed to overcome these issues thanks to the design 
of instrumented particles [13-16]. This was initially developed by our group to study Lagrangian particles 
with a temperature sensitive dependance, there used in Rayleigh-Benard convection [ ]. The approach is 
to instrument a neutrally buoyant particle in such a way that it measures the temperature as it is entrained 
by the flow, and to transmit the data via a radio frequency link to the lab operator. This way, one gains 
access to trajectories for as long as the particle's battery lifetime. In the work reported here, we built upon 
this approach to instrument the particle with a 3D accelerometer such that one gets the accelerations - 
i.e. the forces - acting on a spherical particle in real time and for long trajectories. The instrumented 
particle has been previously tested, benchmarked and validated with the optical technique of Zimmermann 
et al. [ ] , showing a good agreement between the two different measurements of the acceleration. In the 
present work we establish methods to extract physical characteristics of the investigated flows from the 
particle's acceleration signal. 

One further motivation is to gain insights into a flow when direct imaging is not possible, e.g. when deal- 
ing with opaque vessels, non-transparent fluids or granular media. These constraints occur especially in 
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FIG. 1: a) Picture of the instrumented particle (so-called smartPART from smartlNST S.A.S.); b) Possible 
mass distributions of the particle; its inertia consists mainly of a disk and a spherical shell, with different den- 
sity adjustment and imbalance settings by adding Tungsten paste (in blue); experiments are best done with a 
symmetrical repartition of the masses. 

industry where additional bio-medical or environmental constraints arise (the injection of tracer particles 
might be unsuitable and thus prevent any visualization technique). As mentioned above solid particles 
are found to rotate when advected in a highly turbulent flow [12]. We show here that it is possible to 
build quantities depending or not on the particle's rotation and we conclude about flow parameters that 
are directly accessible without any optical measurement. 

The article is organized as follows: first, we present the experimental setup, as well as a brief reminder 
of the technical characteristics of the instrumented particle and the forces it measures (section I). Then, 
we present the new signal processing methods (section II). Finally, we discuss and conclude on this new 
measurement technique (section III). 

I. EXPERIMENTAL SETUP 
A. Instrumented particle 

The device described in the following is designed and built by smartlNST S.A.S., a spin-off from CNRS 
and the ENS de Lyon. It consists of an instrumented particle (the so-called smartPART ®), a spherical 
particle which embarks an autonomous circuit with 3D-acceleration sensor, a coin cell and a wireless 
transmission system, and a data acquisition center (the so-called smartCENTER ®) which receives, 
decodes, processes and stores the signals from the smartPART (see Fig. 1 and Fig. 2). The smartPART 
measures the three dimensional acceleration vector agp acting on the particle in the flow. It is in good 
agreement with other techniques, details can be found in Ref. [18]. The accelerometer consists of a micro- 
electro-mechanical system giving the three components of the acceleration (each of the three decoupled 
axes returns a voltage proportional to the force acting on a small mass-load suspended by micro-fabricated 
springs). From this construction arises a permanent measurement of the gravitational force/acceleration 
g = 9.8 m/s 2 • e z = g ■ e z . Each axis has a typical full-scale range of ±3.6 g = 35 m/s 2 . The sensor has to 
be calibrated in order to compute the physical accelerations from the voltages of the accelerometer. The 
detailed procedure is described in Ref. [ : ]. Concerning the resolution of the smartPART, the uncertainty 
on the acceleration norm is |<x| = Vz^i °f = 0.008 g, with an average noise Oi < 0.005 g on each axis. 

The particle rotates freely and in an a priori unknown way as it is advected by the flow. The instanta- 
neous orientation of the particle can be described by an absolute orientation with respect to a reference 
coordinate system, R(£) [19]. For readability we only write the time reference when necessary {e.g. when 
two different times are involved in an equation) and drop it otherwise. Using this rotation matrix, it is 
possible to express the contributions to the force acting on the particle and measured by the acceleration 
sensor in the lab frame or in the particle frame. The following contributions to the particle's acceleration 
signal agp can be identified: 

(i) Gravity: By construction, gravity g is always contributing to asp- Since the particle is a priori 

oriented arbitrarily in space, g is projected onto all 3 axes. 

(ii) Translation: The forces acting on the particle are projected as the Lagrangian acceleration 

itrans = 'T~o x v : ) onto the sensor. 
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(iii) Rotation: The particle rotates freely around its geometrical center with an angular velocity u>. If 
the sensor is placed by r outside the geometrical center of the sphere one observes the centrifugal 

force: a c f = u> x (uj x r) + ^— u^j x r. According to the technical drawings it is r w 3 mm - e z . 

Experiments on the rotation of the smartPART in a von Karman flow created by two counter- 
rotating impellers show that the angular velocity w of the particle is of the order of the impeller 
frequency /i mp [12, 18]. The rotational forces are of order ruJ 2 ~ r- (in 2 f? mp ) < 0.1 |atrans| and 
have consequently negligible effect. A more detailed analysis measured that the ratio between the 
contribution due to the rotation and the total acceleration is |a c f|/|a tra ns + a c f | < 0.1 [18]. The 
contribution due to the rotation is thus neglected. It has to be noted that by construction of the 
accelerometer and because the circuit is fixed within the sphere, there is no contribution of the 
Coriolis force. 

(iv) Noise &; spikes: In ideal situations the smartPART has a noise of less than 0.005 g for each axis, 
which can be handled by a moving average. Wrong detections appear as strong deviations from 
the signal and are hard to distinguish from high acceleration events due to the turbulent flow or 
contacts with e.g. the impellers. Experiments in different configurations prove the remaining noise 
to be negligible [18]. 



Combining the different terms, and neglecting possible noise and the rotational bias yields: 



a SP w R 



d 2 ' 
9+ W X 



E [9 + "trans] ■ (!) 



The contributions due to gravity and translation are thus entangled by the continuously changing 
orientation of the particle. Since gravity is of little interest here, one has to investigate how common 
quantities such as the mean and the variance (or rms) of the acceleration time series as well as auto 
correlation functions can give information about the particle motion. 

Considering robustness, the smartPART is able to continuously transmit data for a few days. During 
various experiments in a von Karman flow neither contacts with the wall nor with the sharp edged blades 
of the fast rotating impellers damaged its function or shell. Furthermore, the sensor has among other 
things been chosen for its weak temperature dependance; in order to achieve optimal precision of the 
measurements, we calibrate the particle at experiment temperature shortly before the actual experiment. 

Finally, by adding Tungsten paste to the inside of the smartPART the weight of the particle can be 
adjusted such that the particle is neutrally buoyant in de-ionized water at 20 °C. It should be noted 
that the mass distribution inside the particle is neither homogeneous nor isotropic. The particle's inertia 
is best described by a heavy disk of 20 mm diameter (the battery), a spherical shell and patches of 
Tungsten paste. The paste must, therefore, be added carefully to minimize the imbalance of the particle 
(see Fig. lb); otherwise the resulting out-of-balance particle (i.e. with the center of mass not coinciding 
with the geometrical center) induces a strong preferential orientation and wobbles similar to a kicked 
physical pendulum. For a well balanced particle, which rotates easily in the flow, one of the eigen-axes of 
inertia then coincides (approximately) with the z— axis of the accelerometer. The other two are within 
the x — y plane due to rotational symmetry. 



B. Von Karman swirling flow 



We investigate the motion of the instrumented particle in a fully turbulent flow. Namely, we use a 
von Karman swirling flow; in contrast to Ref. [12] the apparatus is here filled with water and develops 
higher turbulence rates. A swirling flow is created in a square tank by two opposing counter-rotating 
impellers of radius R = 9.5 cm fitted with straight blades 1 cm in height (sec Fig. 2). The flow domain in 
between the impellers has characteristic length H — 20 cm = 2R. Blades on the impellers work similar 
to a centrifugal pump and add a poloidal circulation at each impeller. For counter-rotating impellers, 
this type of flow is known to exhibit fully developed turbulence [20] . Within a small region in the center 
the mean flow is little and the local characteristics approximate homogeneous turbulence. However, at 
large scales it is known to be anisotropic [21, 22]. Key parameters of the turbulence at different impeller 
speeds are given in Table I. The two impellers can also be driven co-rotating, creating a highly-turbulent 
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FIG. 2: a) Picture of the von Karman swirling flow with the instrumented particle inside; b) Sketch of the ex- 
perimental setup with the apparatus and the smartPART transmitting acceleration signals to the smartCENTER; 
c) Sketch of the global structures that can be found with the two co- and counter-rotating regimes. 
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TABLE I: Key-parameters of the counter-rotating flow configuration. The integral time scale is defined as Tint = 
l//im P and the integral length scale is estimated to be L[ nt = 3 cm. We use the following definition for the 

Reynolds numbers: Re = 2ttR 2 fim P /u and R\ ~ \J ^ ■ 27rL? nt /j mp . Note that the particle explores the whole 
apparatus, where the flow is known to be inhomogeneous and anisotropic. Thus, R\ and the Kolmogorov scales 
are only rough estimates. For comparison: co-rotating impellers yield an energy injection rate which is half of the 
energy injection rate of counter-rotating impellers at the same impeller frequency. 

flow inside the vessel with one pronounced persistent global vortex along the axis of rotation. Close to 
the axis of rotation the mean flow is weak, followed by a strong toroidal component and an additional 
poloidal circulation induced by blades on the impellers. The energy injection rate is a factor of 2 smaller 
than for counter-rotating impellers. This means that at the same impeller frequency, /imp, co-rotating 
driving creates less turbulence than counter-rotation, but the flow is still highly turbulent [23] . Although 
the vortical structures near the disks are comparable (see Fig. 2c), the co- or counter-rotating regimes 
yield well distinct global structures in the center of the vessel. The two regimes are used to compare the 
signals obtained by the instrumented particle in two very different flow configurations. In addition, the 
co-rotating serves as a test case for persistent, large vortex structures as they are found in mixers with 
only one impeller. 

II. ACCELERATION SIGNALS 

Fig. 3 shows two sample time-series of the three components of the acceleration measured by the 
instrumented particle in the von Karman flow, superimposed with the norm of the acceleration. Two 
different frequencies of the impellers are presented here: 1Hz and 4 Hz. The mean value of the norm 
fluctuates around g = 9.8 m/s 2 , indicating that gravity is always measured by the accelerometer. 
Furthermore, the fluctuations of the norm increase with the impeller frequency. It is, however, 
difficult to compare the three components of the acceleration, either between each other or for dif- 
ferent impeller frequencies. This is mainly due to the measurement of the gravity that is randomly 
projected on the three axes of the accelerometer as the particle rotates in the flow. It results in 
signals containing both contributions from the gravity and the particle's translation, with no straight- 
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forward method to separate them. In other words, contrary to other methods (e.g. particle tracking 
velocimetry) it is not possible to obtain directly the characteristics of the particle motion. Hence, one 
needs to post-process the data to derive information about the statistics and the dynamics of the particle. 



A. Analysis of the raw signal asp 



Fig. 4 shows different results of a basic statistical analysis of the acceleration signals, namely the PDFs 
of the three components and the norm of the acceleration for different impeller frequencies, and the 
fluctuating and mean values of the acceleration as a function of the impeller frequency. The accelerometer 
used in the smartPART saturates if one of the acceleration components exceeds ±3.6(7, we exclude these 
points from the analysis. This removal diminishes the observed acceleration and the bias increases with 
the forcing. In the case of Fig. 4, almost 3% of all data points were removed at /; mp = 4 Hz, which is two 
orders of magnitude higher than for /j mp = 1 Hz. Looking at the PDFs of the acceleration for a given 
impeller frequency (Fig. 4a), one can see that the three components give similar results for a wide range 
of acceleration values. However, the PDFs are very different from one frequency to another. Whereas 
at low impeller frequencies the PDFs are skewed and shifted, they become centered and symmetric 
with increasing impeller frequency. This evolution in shape can be explained by the particle's mass 
distribution and imbalance. Although the particle is carefully prepared, its moment of inertia is not that 
of a solid sphere and the particle's center of mass does not perfectly coincide with its geometrical center. 
Consequently, the particle becomes slightly out-of-balance, with a preferred orientation at low impeller 
frequency: the peaks then correspond to the projection of g on the axes in this preferred orientation and 
fluctuations around it. When the impeller frequency (and consequently the turbulence level) increases, the 
particle is able to explore all the possible orientations, meaning g is randomly projected in all directions, 
and the asymmetry disappears. The PDFs of the norm |ogp| (Fig. 4b) also show this difference in shape, 
with a clear peak near the value g (again, gravity is always measured by the 3D accelerometer), but with 
a narrow strong peak at low impeller speed and a more stretched PDF at high impeller speed. 

The evolution of the fluctuations of acceleration (rms asp) with the impeller frequency is given in 
Fig. 4c. Only one component of the acceleration is presented here for readability, since no preferred 
direction in any of the axes was found. This results in all three components rms values having the same 
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FIG. 4: Top: Probability density functions of a) the components and b) the norm of the acceleration asp for two 
different impeller frequencies. For readability, the PDFs of acceleration components at 4 Hz have been arbitrarily 
shifted. Bottom: Evolution of c) rmsasp and d) (asp) with the impeller frequency /i mp ; filled symbols (o) 
indicate co-rotating impellers. In all cases the particle explored the flow for a sufficient amount of time for the 
statistics to converge. In good agreement with Eq. (2), |{asp)| continuously decreases from lg to Og as the impeller 
frequency increases. 



behavior and amplitudes. One can see that surprisingly, the fluctuations of acceleration increase linearly 
with the frequency. That is in contrast to dimensional arguments that tell ot^ rans oc 

/knp- Moreover, it is 

not possible to distinguish between the co- and counter-rotating regimes of the impellers. 

Fig. 4d shows (asp) as a function of the impeller frequency, /imp, and the forcing. As expected the 
mean accelerations are becoming smaller with increasing impeller frequency. Indeed taking the average 
of Eq. (1) yields: 

(a SP ) = (Rg) + (Eatrans) • (2) 

If the particle explores continuously all the possible orientations, the mean vanishes; whereas a fixed 
orientation (i.e. no rotation) yields (asp) = E(9 + (okrans))- This is what is observed in Fig. 4. In 
the case of weak turbulence (i.e. for smaller values of /imp), the mean acceleration gives an estimate 
of gravity: (asp) ~ (Es)- However, for stronger turbulence and even if the mass distribution slightly 
induces a preferred direction, the particle can rotate freely around this axis, resulting in a vanishing mean 
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acceleration when the impeller frequency increases: (agp) ~* 0. In the latter case, contacts with impellers, 
walls, eddies, etc also help overpowering any preferred direction easily and force the particle to rotate. It 
can be noted that again, it is not possible to distinguish between the co- and counter-rotating regimes. 
Furthermore, the variance of a component agp,i of agp depends strongly on its mean value, (asp,i)- As 
explained before, gravity renders (<xsp,i) non- negligible. Additionally, we observe for weak turbulence 
levels (/imp ;$ 1 Hz) that particles are able to stay in an orientation for several seconds. Hence, a global 
mean of the complete time-series is not a meaningful quantity. 

The direct study of the raw acceleration signal, agp, only allows to conclude whether the particle rotates 
or not. It does not permit to disentangle the contributions of the gravity and the particle translation, 
and subsequently to have any precise insight on the flow. Other methods adapted to this problem are 
thus needed to extract informations from the instrumented particle related to its motion. 



B. Moments of the acceleration due to the particle's translation 

In confined flows and provided the statistics are converged, it is (a tra ns) = 0. One is, therefore, 
interested in the PDF of atrans- Although, we mentioned we don't have direct access to atrans an d its 
PDF, we can compute the even (central) moments of its PDF. 

The variance of agp is 

(Ogp) = (lIM 2 ) + (|]|a trans | 2 ^ +2(ES'-|| atrans) 
= 5 2 + (<4ans)+2g(a z ), 



(3) 



where a z = e z ■ a tra ns- It should be kept in mind, that each axis of the smartPART's accelerometer is 
limited to ±3.6 g, and possible events of higher acceleration are therefore not included in the analysis. 
The PDF of |agp| 2 for different impeller frequencies is shown in Fig. 5a. As expected, a peak is 
clearly observed at g 2 . One can also see that there are breaks in the slope at |agp| 2 ~ (3.6 g) 2 and 
| agp | ~ 2(3.6g) 2 , corresponding to cases where one or two axes would saturate. Some information is 
inevitably lost, and to investigate the behavior at large /i mp , the sensor would have to be replaced with 
a different model supporting higher accelerations. 

If the particle is neutrally buoyant and the flow is confined, one expects (a z ) — 0. We therefore obtain 
an estimate of the standard deviation (rms) of a trans : 



V^ans) = V / (4p)-5 2 - ( 4 ) 



a rms is independent of how gravity is projected on the axes of the accelerometer (in other words 
it is insensitive to the particle's absolute orientation). Nevertheless, a bad calibration (e.g. caused 
by longterm drift or strong temperature change) can introduce a systematic offset to a rms . Nev- 
ertheless, this bias can be minimized by calibrating the thermalized smartPART before the actual 
experiment. Fig. 5b depicts the evolution of a rms with the impeller frequency. In agreement with 
dimensional analysis, a rms (/i mp ) describes a parabola, although one could also argue that a rms seems 
linear with /i mp for /j mp > 3 Hz. As illustrated in Fig. 5a, this is caused by saturation of the ac- 
celerometer, which cuts off/underestimates high acceleration events present at these high turbulence level. 

Similarly to the variance, one can estimate the fourth central moment of a trans . It is 
< I^spI 4 ) - ( [g 2 + a t 2 rans + 2 g a z ] [g 2 + a t 2 rans + 2 g a z ] ) 

(5) 

= g 4 + { Kansl 4 ) + 2.g 2 <a 2 . ans ) + Ag 2 (a 2 z ) + Ag 3 (a z ) + Ag (|a trans | 2 a,) . 

Assuming no preferred direction in a t rans, as found for small particles in a windtunnel [4] and verified 
for solid particles of size comparable to the integral length scale in the same apparatus [ ], one has 
Ag 2 (a 2 ) ~ 4/3 g 2 a 2 rans . Again, the terms (a z ), A g 3 (a z ) and Agl |a trans | 2 a z \ are expected to vanish in 
the case of confined flows. Eq. (5) then simplifies to 

<|asp| 4 )«g 4 + <K an s| 4 ) + f 5 2 «r 2 m s- (6) 
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0123401234 
impeller frequency [Hz] fi mp [Hz] 

FIG. 5: Moments of a tra ns- a) PDF (a|p) at different impeller frequencies. The 3 vertical lines mark gravity and 
the saturation of one or two accelerometer axes, b) and c) : RMS and fourth moment of atrans as a function of 

the impeller frequency. a rms (/i mp ) = J (ja S p(/i mp )| 2 ) - g 2 . 



The flatness, -^(atrans): is defined as: 

s (|atra„s| 4 ) ( |a SP | 4 ) - .g 4 - f ff 2 ^ 

F (Otraiu.) = — ~ 2 72 = ~i " • ( 7 ) 

\ a trans/ rms 

As shown in Fig. 5c we observe a flatness of the order of 10 in our von Karman flow, which is 
close to the flatness obtained in the case of much smaller particles [6] and to our finding for solid 
particles of similar size [ ]. The uncertainty in the flatness can partially be attributed to an un- 
certainty in g and stems from the resolution, noise and measurement range of the smartPART, but 
also from the particle's weak drift. It is furthermore biased by contacts with impellers and walls. 
More surprisingly, the flatness decreases with the forcing. This decline is again due to the limited 
measurement range of the accelerometer used: at high accelerations the sensor saturates and thereby 
sets PDF (atrans | |osp,i| > 3.6<?) = 0. Since the flatness is the fourth moment of the PDF and as such 
highly sensitive to strong accelerations, we find a decrease whereas solid large spheres in the same flow 
have an increasing flatness [ ] . We also conclude that calculating moments of higher order is out of reach. 

It is remarkable, that based only on the second and fourth moment of atrans one cannot clearly distin- 
guish between a counter-rotating and a co-rotating flow although these two forcings induce two clearly 
different large scale flow structures. Similar behavior has been found for solid spheres of comparable size 
in the same flow [ I I], too. It should be noted that the energy injection rates for the two ways of driving 
the flow differ by only a factor of 2: the co-rotating forcing is highly turbulent, too. In addition, in 
vicinity of the disks the flow has a strong contribution of the centrifugal pumping of the blades on the 
impellers and the flow configurations are comparable in that region. 
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C. Auto-correlation functions 



In order to distinguish between the two regimes, we now turn to the auto-correlation of the acceleration 
time-series to estimate correlation time scales of the flow. Ideally, one would want to compute the auto- 
correlation of the translational force, e.g. (atrans (t) ' atrans (t + t)) , but again the constantly changing 
orientation of the smartPART blocks any direct access to atrans (t) and quantities derived thereof. We 
therefore need to find quantities, which are either not altered by the orientation of the smartPART or 
extract information on its rotation. 



1. An Auto-correlation invariant to the rotation of the sensor 



+ 




10 



10 10 
correlation time T [s] 
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FIG. 6: Top: Rotation-invariant auto-correlation function C m0 v (t) (Eq. (8) after rescaling); the auto-correlation 
significantly differs between counter- and co-rotating impellers. Bottom: Rotation-sensitive auto-correlation func- 
tion (asp(i) ■ asp(i + r)) / g 2 . In all cases the impeller frequency is 1.5 Hz. A logarithmic scale is chosen for the 
abscissae to display both short and long time contributions to the correlations. 



In the spirit of Eq. (3) and Eq. (5) one can construct the auto-correlation function of the magnitude 
of agp. It is: 

(|a S p(i)| 2 |a S p(t + T)| 2 ) 

= < [.9 2 + a 2 tlans (t) + 2g- a trans (i)] [g 2 + a 2 . ans (i + r) + 2 g • a trans (i + r)] ) 

= (|at r ansW| 2 |atrans(t + T)| 2 ) + ff 4 + 2 5 2 ( Q 2 rans ) + 4g 2 (a z (t) a z (t + r)) (8) 

+ 4g 3 (a z ) + 2g (\a trans (t + t)| 2 a z (t)) + 2g <J|a trans (i)| 2 a z (t + T )} 

w (|atra„s(t)| 2 |a trans (i + t)| 2 ^ + g 4 + 2 g 2 a 2 . nls + 4g 2 (a z (t) a z (t + r)) . 

Again, the terms containing a z = e z ■ atrans a re expected to have zero mean. However, the last term on 
the right-hand side of Eq. (8) does not vanish for r pa 0, becoming Ag 2 (a z (t) a z (t + t)) = Ag 2 (^\a z \ 2 ^j. 

Assuming no preferred direction in a trans this can be approximated as 4/3 g 2 a 2 nls . In contrast to Eq. (8), 
we preferably compute the autocorrelation of the fluctuations around the mean (a|p). Hence, the 
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autocorrelation of the norm can be negative. We further normalize the auto-correlation such that it is 1 
at t = 0. After rescaling, we refer to this rotation- invariant quantity as C mov (r). 

Fig. 6 displays C mov (r) for the co- and counter-rotating regime at an impeller frequency of 1.5 Hz. The 
auto-correlation of the counter-rotating forcing is well approximated by a sum of exponential decays. In 
contrast thereto, we observe that co-rotating impellers correspond to an auto-correlation function showing 
a damped oscillation, i.e. the smartPART observes the longer coherence in the large scale motion of 
the flow. This is in agreement with Eulerian measurements, where pressure probes were mounted in a 
von Karman flow: whereas the counter-rotating flow produces typical pressure spectra, the same probe 
in the co-rotating case yields a spectrum which peaks at multiples of the impeller frequency. Similar 
behavior has been reported for the magnetic field in a von Karman flow [24] filled with liquid Gallium 
(in that particular case the co-rotating regime was created by rotating only one impeller). 

Summing up, / |agp(i)| 2 |agp(f + r)| 2 ) is insensitive to the particular tumbling/rotational dynamics of 
the particle and it gives necessary information to determine the type of flow. We also checked that this 
result is not altered by a possible imbalance of the particle. 

2. An Auto- correlation related to the tumbling of the particle 

One can further focus on the rotation of the particle by considering the product: 
(a SP (t) ■ a SP (t + t)) 

= ([E(i) (9 + (<))] ■ [R(i + r) (g + ^trans (t + r))]) 
= <[R T (i + r) .g(t) (g + a t 

(*))] • i(9 + & trans (t + r))}) 

^trans 

(*)] (t + r)) + ([T(t,r)g) 

' ^trans (t + T )) + ([T(i, T)a trans (i)] g) 

» 9 2 (e, • [W, r)e z }) + ([Tji, r) • a trans (t)} a trans (t + r)) , 

(9) 

where the term T(t, r) = R T (< + r)R(t) is a rotation matrix related to the instantaneous angular 
velocity w of the particle as explained in Ref . [10] . Again, the two terms containing products of g and 
a vanish if the particle is neutrally buoyant. The term g 2 (e z ■ [T(£,r)e z ]) is related to the tumbling of 
a spherical particle [ ]. In contrast to the other auto-correlation Eq. (8), one cannot subtract a mean 
value prior to computing (asp(t) ■ asp(t + r)) . To estimate the ratio between g 2 (e z ■ [T(t, r)e 2 ]) and 
(LT(* ; 7")atrans(*)] ' atrans(i + t)) it is helpful to normalize by g 2 . 

If ([T(t, t) ■ a t rans(i)] a trans(^ + 7")) becomes uncorrelated, (agp(t) ■ asp(t + r)) does not necessarily van- 
ish. If uncorrelated (r 3> Ti nt ): 

(a SP (t) ■ a SP (t + r)> = g 2 (R (t) e z ■ E (t + r) e z ) 

=g 2 (m*) 2 >° 111)1 

That means (asp(t) • asp{t + r )) approaches a plateau whose height is determined by the average orienta- 
tion of the particle. In analogy to C mov , one can then subtract g 2 (Re z ) and rescale (asp(t) • o,sp(t + T )), 
which is termed C tum b (t) in the following. 

The lower plot in Fig. 6 depicts (asp(t) • asp(t + r)) /g 2 for the two forcing regimes at an impeller 
frequency of 1.5 Hz. In agreement with Eq. (10), a plateau is reached for both forcings. To investigate 
the role of the plateau we plot the auto-correlation of the particle for increasing /j mp in Fig. 7. For 
/imp £ 1 Hz one finds little change with the plateau at almost 1. For /; mp w 2 Hz the plateau drops but 
is still non-zero. The value of the plateau diminishes with further increase in /i mp . At the same frequency 
range we observe that the PDFs of the components of agp become centered and symmetric (cf. Fig. 4). 
Ctumb (t) - i.e. Eq. (9) after subtracting the plateau and rescaling - accesses the fluctuations around a 
mean value and is shown in the bottom plot of Fig. 7. We find that C tum b (t) evolves from a long-time 
correlated oscillatory shape at low impeller speeds to an exponential decay at high /i mp . 
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FIG. 7: Top: Dependence of (asp(t) ■ asp(t + r)} / g 2 on the impeller frequency for counter-rotating impellers; 
Bottom: Same data after subtracting the plateau and rescaling. 



3. Time scales 



The autocorrelation functions contain time-scales which are related to the movement of the particle in 
the flow. We identify two scenarios for the autocorrelations (C mov (t) and C tU mb (t)): as illustrated in 
Fig. 6, they are either conducting a weakly- or a critically-damped oscillation. With increasing turbulence 
level the oscillation is gradually changing towards the critically damped case and for high propeller speeds 
(/imp > 2.5 Hz) no damped oscillation is observed (cf. Fig. 7). In order to extract meaningful time-scales 
we, therefore, fit two test-functions to each autocorrelation function. The functions are the transient 
solution of a weakly damped harmonic oscillator: 

f w (r) = a exp (-r/r corr ) • sin (27r/ osc + (f> ) , (11) 

and of a critically damped one: 

fd(r) = exp (-t/t co1t ) ■ (a + ai r) . (12) 

Tcorr,/osc and ao,ai,<po are fit-parameters. We return T cor r and (if available) / osc from the test-function 
which performs better in approximating the autocorrelation. C mov (r) and Ct um b (t) access motion and 
tumbling of the particle, respectively, and yield thus different time scales. For the oscillatory case, / osc 
contributes additional details on the particle's motion. 

Fig. 8 shows T corr as a function of the impeller frequency and driving. We find that both rotation- 
invariant (Eq. (8)) and rotating-sensitive (Eq. (9)) function find very similar correlation times in the 

3/2 

counter-rotating configuration. Moreover, r corr of the particle follows roughly a / imp power-law as 
suggested by the scaling of the Kolmogorov time scale (it is t v cx e" 1 / 2 and e oc f? mp ). Furthermore, r corr 
obtained from the rotation-sensitive function is independent of the way we drive the flow. In contrast 
thereto, the rotation-invariant function gives correlation times, r corr , which are larger and only little 
dependent on the impeller frequency if the impellers are co-rotating. That means one can distinguish co- 
from counter-rotating forcing by comparing the timescales of the two auto-correlation functions. 

Concerning the oscillation frequency / osc (not shown in figure), we find that the rotation-sensitive 
autocorrelation senses the tumbling/wobbling of the particle, which is directly related to the particle's 
imbalance and independent of the flow. The rotation-invariant autocorrelation on the other hand shows 
an oscillation frequency following the impeller frequency with / osc ~ §/imp- 
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1 2 3 

impeller speed fi mp [Hz] 



FIG. 8: Characteristic flow time scale r corr determined from rotation-invariant, C mov (t), and rotation-sensitive, 
Ctumb (t), auto-correlation function. Whereas both autocorrelation functions yield comparable T cor r for counter- 
rotating driving, we find that their results are well distinct in the case of co-rotating impellers. The dashed line 
indicates a f^/ 2 power-law, as suggested by the scaling of the Kolmogorov time scale (r v oc e~ 1//2 and e oc /^p). 



III. DISCUSSION 

After briefly presenting the working principle of an instrumented particle measuring Lagrangian ac- 
celerations, we established a mathematical framework based on statistical moments and auto-correlation 
functions to analyze turbulent flows from the particle's signals. In particular, we developed methods 
which are either invariant or sensitive to the rotation of the particle and its sensor in the flow. These 
methods perform well within the wide range of tested turbulence levels. With a smartPART one gets 
access to correlation time scales of the flow, as well as the variance and flatness of the (translational) 
acceleration. Comparing the rotation-sensitive and the rotation-invariant autocorrelation allows distin- 
guishing between different flow regimes, notably detecting long-time correlated large vortex structures as 
shown here with the co-rotating forcing of a von Karman flow. In contrast to particle tracking methods 
the instrument particle returns one long trajectory instead of many short realizations. To that extent, it 
has to be noted that we limited our analysis to the extraction of global flow features. In order to follow 
the evolution of a slowly changing flow in time, these methods can, however, be extended to sliding win- 
dows. Work on adaptive filtering techniques is ongoing, in particular we are testing the Empirical Mode 
Decomposition, which might be able to separate the different contributions of the signal and thereby get 
even deeper insight into the flow. 

We emphasize that after usage the particle can be easily extracted from the flow and then be reused and 
that by virtue of the developed mathematical framework no optical access is needed. The instrumented 
particle can therefore shed some light into flows that are not or hardly accessible up to now. Due to its 
continuous transmission one flow configuration can be characterized within ~ 30min. This technique is 
an interesting tool for a fast quantification of a wide range of flows as they are found both in research 
labs and industry. 
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